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Abstract. We present a detailed study on the kinematics of metal-rich stars with and without planets, and their relation with 
the Hyades, Sirius and Hercules dynamical streams in the solar neighbourhood. Accurate kinematics have been derived for all 
the stars belonging to the CORALIE planet search survey. We used precise radial velocity measurements and CCF parameters 
from the CORALIE database, and parallaxes, photometry and proper motions from the HIPPARCOS and Tycho-2 catalogues. 
The location of stars with planets in the thin or thick discs has been analysed using both kinematic and chemical constraints. 
We compare the kinematic behaviour of known planet-host stars to the remaining targets belonging to the volume-limited 
sample, in particular to its metal-rich population. The high average metallicity of the Hyades stream is confirmed. The planet- 
host targets show a kinematic behaviour similar to that of the metal-rich comparison subsample, rather than to that of the 
comparison sample as a whole, thus supporting a primordial origin for the metal excess observed in stars with known planetary 
companions. According to the scenarios proposed as an explanation for the dynamical streams, systems with giant planets could 
have formed more easily in metal-rich inner Galactic regions and then been brought into the solar neighbourhood by dynamical 
streams. 
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1. Introduction 

Precise spectroscopic studies have revealed that stars with 
planets seem to be particula rly metal-rich when comp ared 
I with "single" fi eld dwarfs (e.g. Gonzalez 1997; Gon zalez et alJ 
l200lUSantos et alJ200lll2004bH2005i) . This represents the first 
strong link between the planetary companion and its parent 
star. In the last few years we have published a series of pa- 
pers on chemical abundances in stars with known giant plan- 
ets. Abundances of volatile and refractory elements in a large 
set of planet-host stars an d in a volume-l i mited compariso n 
sample were presented in |E cuvillon et al. (2004a b, 2006b); 
iBodaghee etallfeOoMBeirao et alJJ2005l) : lGilli et al) fe006). 
The volatile/ refractory abundance ra tios have recently been in- 
vestigated in lEcuvill on et alj J2006al). Be and Li light elements 
were also analysed JSantos et alJl2002al l2004al [israelian et alJ 
120041) . We are now interested in investigating the kinematic be- 
haviour of planet-host stars. 

The kinematics of stars in the solar neighbourhood can 
give fundamental information for our understanding of the 



structure and evolution of the Milky Way. ESA's astrometric 
satellite Hipparcos (ESA 1997) has provided us with accu- 
rate positions and trigonometric parallaxes, as well as abso- 
lute proper motions, for a large and homogeneous sample of 
tens of thousands of stars near the Sun. This has offered the 
opportunity to investigate the velocity distribution in the solar 
neighbourhood, not only for early-type stars, but also for the 
old population of the Galactic disc. S everal studies have been 
performed on this topic since (e.g.. iDehnen & Binnevl 19981: 
| Dehnenlfl993 lAlcobe & Cubarsil l2005t IChereul et alJ " 
Hoogei 



|1999| 




Asiain et al,| 


1999 
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The kinematic distribution of stars in the Solar neighbour- 
hood is far from smooth. The inhomogeneities can be spatially 
confined groups of young stars (e.g. young clu sters), or spa- 
tially extended groups with the same kinematics. Eggen lll994l) 
defined a "supercluster" (SC) as a group of gravitationally un- 
bound stars that share the same kinematics and may occupy ex- 
tended regions in the Galaxy, and a "moving group" (MG) as 
the part of the supercluster that enters the solar neighbourhood 
and can be observed all over the sky. 
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According to the standard interpretation, stars in a mov- 
ing group were formed simultaneously in a small phase-space 
volume, so that we can still observe a stream of young stars 
with similar velocities. On the other hand, old stellar popula- 
tions should be completely mixed and present a smooth distri- 
bution function, due to phase mixing and scattering processes. 
Actually, there are basically two factors acting against the per- 
sistence of a moving group in the general stellar background: 
Galactic differential rotation, which tends to spread them out 
very q uickly in the direction of the Galac tic rotation dWoollevI 
1960), and disc heating (e.g. lLace vlll99ll) . which prod uces an 
increa se in the velocity dispersion of disc stars with age Wielenl 
ill 9771) reported that a typical star in the solar neighborhood 
would forget its initial peculiar velocity within 2 x 10 8 years. 

However, several works have obtained results that raise 
problems with this interpretation. On the one hand, it has been 
reported that early-type stars belonging to the same superclus- 
ter span a wide range of ag es, in co ntradiction with the assump- 
tion of a common origin (Chere ul et alJ ll998. 1999). On the 
other hand, it is striking to verify that some of the classical 
MGs are several 10 8 yr old, and that significant dum piness in 
veloc ity space has been reported for late-type stars iPehnenl 
Il998h . Other mechanisms of different nature have been pro- 
posed as responsible for the substructures observed in veloc- 
ity space. Streams could be the remnants of merger events be- 
tween the Mi lky Way and a satellite galaxy, as recently pro- 
posed by Nav arro et alJj2004l) for the Arcturus group, although 
the chance that two of them have left such impo rtant signature 
in the disc near th e Sun is statistically very low(Famaev et alJ 
2005a). Moreover, Rab oud et al 1 (1 19981) reported that the local 
anomaly found in the (u, v) plane is caused by intermediate- 
to high-metallicity stars, while a merger would involve rather 
metal-poor stars (Dehnen 2000J). 

Another possibility is that purely dynamical mechanisms, 
such as non-axysimmetric components of the gravitational po- 
tential (for example, the rotating galactic bar or spiral waves), 
are the cause of substructures in the velocity distribution (e.g . 
iDehnenll 19981 iRaboud etal.ll 19981 ISellwood & Binnevl Eo02). 
Stars in an MG could have formed in the inne r part of the disc 
and have become trapped in resonant orbits. Dehnen (1999) 
associates the bimodality of the local velocity distribution of 
late-type disc stars, lately identified as the Hercules stream, 
with high-meta llicity stars thrown out from the inner disc, and 
Dehne 1 ll2000l) concluded that outer Lindblad resonant scat- 
tering off the Galactic bar is presently the only viable expla- 
nation for it. Other simulations have tested the effect of other 
non-axysimmetric components, in particular transient spiral 
waves, in the Galacti c disc (e.g. ISellwood & Binnevl f2002: 
iDe Simone et al.l2"0 04). with the result that they may cause ra- 
dial migration near their corotation radius. 

It is not unlikely that known exoplanetary systems could 
also have suffered such radial displacements. The typical con- 
tributors to the kinematic str uctures in the solar neighbour- 
hood are metal-rich old stars jRaboud et alJll998i) . The well- 
known metal-rich nature o f planet-host stars ( e.g. iGonzaled 
1 19971 iGonzalez et all 1200 it ISantos et aljEoOll l2004bl 120051) 
might thus be related with dynamical streams. It is probable 
that planet-host stars may have formed in protoplanetary clouds 



with a high metal content, such as those located in inner re- 
gions of the Milky Way, and then be brough t closer to the 
Sun by dynamical streams dFamaev et alJ2005bl) . lFamaev et alJ 
d2005b|) have found an in dicative clump of stars hosting planets 
from lSantos et al 1 d2003l) in the Hyades streams. Several works 
have investigated t he kinematics of stars with extrasolar plan- 
ets re.g.lGonzalezlll999t lReidll2002l barbieri & Grattonll2 002: 



•• - ^^^gZ]^^^ _ .... . ... . — • — - — ■ 

ISantos et alJl2003f) . without finding any significative peculi 
ity. 

Two main scenarios have been proposed to explain the 
metal excess found in stars with exoplanets. The "self- 
enrichment" hypothesis suggests that the infall of a large 
amount of protoplanetary rocky material on to t he parent 
star may have enhanc ed its primo rdial metallic i ty dGpnzalezl 
Il997l iTaughlinlEoOoT) . In contrast, ISantos et alJ d2000l 12001 1) 
have proposed a "primordial" origin for the iron over- 
abundance of planet-host stars, so that planetary system 
formation would be triggered by a high metal content 



suits supporting the "primordial" 


scenario (Pinsonneault eLajJ 


200 it ISantos et alJl2004bt iReidl 


2002; Sadakanget al.ll2002l 


Ecuvillon et alJl2004albl l2006blal 


iGilli et al.ll2006b. However. 



some clear evidence of pollution events have been re ported 
dGonzaledll998l Laws & Gonzalez! EobH ISmifh et alJbOQll 
Israelian et al. 2001, 2003), so that the "self-enrichment" sce- 
nario cannot be completely discarded dVauclaiJ2004l) . 

This article studies the kinematics of planet-host stars and 
their relation with the dynamical streams in the solar neigh- 
bourhood and investigates whether it corresponds to the be- 
haviour of the metal-rich population or not. The idea is that if 
the metal-rich nature of planet-host stars is intrinsic, then they 
should be kinematically identical to the metal-rich comparison 
sample, while if the metallicity excess is due to planet-related 
enrichment, then they should be kinematically identical to the 
total comparison sample. 

We present accurate kinematics for all the st ars belong- 
ing to the CORALIE planet search survey (see llJdrv et all 
2000). Precise radial velocity measurements and CCF pa- 
rameters were extracted fr om this unique database, whereas 
HIPPAR COS dESAl 1 19971) parallaxes and photometry, and 
Tycho-2 (Hcig et alJl2000 ) proper motions were adopted. The 
location of stars with planets in the thin or thick discs is anal- 
ysed using both kinematic and chemical constraints. We study 
the kinematic behaviour of known planet-host targets with re- 
spect to the remaining stars belonging to the kinematically un- 
biased volume-limited sample and we investigate possible links 
between the presence of planets, the stellar metallicity and the 
dynamical streams Hyades, Sirius and Hercules. 

2. Data and analysis 

We obtained kinematics for all the stars belonging to the 
CORALIE planet search survey. The overall sample consists of 
about 1650 unevolved solar-type stars closer than 50 pc, with 
precise parallaxes (<x n < 5 mas) and spectral types between F8 
and Ml in the HIPPARCOS catalogue. A further selection was 
applied to exclude binaries (cr RV <500ms~'). We used paral- 
lax measurements and photometric data from the HIPPARCOS 
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Fig. 1. Probability distributions of location in the thin disc (up- 
per panel) and the thick disc [central panel), and the thick- 
disc-to-thin-disc relative probability [lower panel). Diamonds 
and crosses represent stars with and without planets, respec- 
tively. In the lower panel, the solid lines indicate the limits for 
reliable membership to each population. 

Table 1. Kinematic parameters adopted for the Hyades, 
Sirius and Hercules velocity ellipsoids. The mean velocities 
[Uo,Vo,Wo) are relative to the LSR. All values are expressed 
in kms -1 . 
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Fig. 2. Thick-disc-to-thin-disc relative probability as a function 
of [Fe/H]. Diamonds and crosses represent stars with and with- 
out planets, respectively. The solid lines indicate the limits for 
the reliable membership to each population. 
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Fig. 3. Upper panel: Metallicity distributions of planet-host 
and comparison stars (solid and dotted lines, respectively). 
The mark indicates the median of the planet-host distribution. 
Lower panel: [Fe/H] vs. stellar mass (in solar mass units) for 
planet-host [diamonds) and comparison [crosses) stars. 

Table 2. Number of targets falling into the velocity ellipsoids 
adopted for the Hyades, Sirius and Hercules streams. 
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catalogue dESAll 19971) . while coordinates and stellar motions 
were extracted from Tycho-2 catalogue facte et alfc OOOl. since 
they improved the HIPPARCOS data by combining several new 
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catalogues. The CORALIE database provided us with radial 
velocity measurements of a unique precision, as well as with 
the average width and depth of the cross correlation function 
(CCF). We derived Galactic space-velocity components U, V 
and W with respect to the Local Standard of Rest (LSR), adopt- 
ing the stan dard solar motion (U^V^W p.) = (10.00, 5.25, 7.17) 
kms 1 from lDehnen & Binnevl(ll998l) . 

2.1. Thin/thick disc 

We selected thin and thick di sc stars from our sam ple following 
the procedure described by Ben sbv et alJ J2003). We assumed 
that the Galactic space velocities of each population (thin disc, 
thick disc and halo stars) have Gaussian distributions. The char- 
acteristic velocity dispersions (cr f/ , cr v , <tw) and asymmetric 
drift (Vasym) were taken from lBensbv et all 1 2003b . as well as 
the observed fractions of each population in the solar neigh- 
bourhood (X). The probability distribution of belonging to a 
specific population (P) is thus given by the Gaussian probabil- 
ity distribution multiplied by the corresponding observed frac- 
tion (X). Stars with P( thick disc)/P(thin disc) > 10 and < 0.1 
will be reliable members of the thick and thin disc, respectively. 
Figure[0shows the probability distributions of belonging to the 
thin disc (upper panel) and to the thick disc (central panel), and 
the thick-disc-to-thin-disc relative probability (lower panel). 
Chemical abundances for CORALIE stars alrea dy spectroscop- 
ically studi ed were extracte d from lEcuvillon et alJ J2004alrJ. 
2006b) and Gilli et al ] J2004 and combined with the kinematic 
information in order to analyse the relation of the trends [X/Fe] 
vs. [Fe/H] with the thin and thick discs. 

All the metallicities were d erived using the cali bration of 
the CCF surface with [Fe/H] bv lSantos etail J2002b). The dis- 
tribution of the thick-disc-to-thin-disc relative probability as a 
function of metallicity is shown in Figure|2] We also computed 
stellar masses by a comparison of the star po sitions in the CM 
diagra m with the calculated isochrones of IPietrinferni et all 
(2004]). A ges for field stars are notoriously difficult to deter- 
mine, and suffer from large biases that make t hem un suitable 
for this kind of study, as discussed bv IPont & Ever! (|2004). 
However, mass can be used as an age surrogate for a first-order 
separation of young and old stars. Masses can be determined 
quite precisely from stellar evolution models and are less sub- 
ject to biases. Figure [5] shows the metallicity distributions for 
planet-host and comparison stars (upper panel), and the plot of 
[Fe/H] vs. mass (lower panel). 
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Fig. 4. Density distribution of Galactic space velocities pro- 
jected on the (V, U) plane, for the whole comparison sample. 
The uncorrected density distribution, the subtracted Gaussian 
distribution and the corrected density distribution are presented 
in the upper, central and lower panels, respectively. The field 
ellipsoid is overplotted with the blue line, while the Hyades, 
Sirius and Hercules ellipsoids are indicated by the red, blue 
green and yellow lines. 



2.2. Analysis of the (U, V, W) space 

The field ellipsoid was defined by computing its inertial tensor: 
the eigenvectors of the inertial tensor would indicate the main 
directions (U',V',W) of the field ellipsoid in the (U, V, W) 
space. We adopted 2 cr of the data along the directions U', V 
and W as the field ellipsoid semiaxes. We then derived the 
density distribution of the Galactic space velocities by the fifth 
nearest-neighbour method. This method assigns to each point 
a den sity value, according to its surroundings (see ISilver man 
1996, section 5.2). 



A field without "structures" is supposed to follow a 
Gaussian distribution. Therefore, this Gaussian "background" 
has to be removed to study the remnant "structures". We thus 
subtracted from the field density distribution the density distri- 
bution of the Gaussian corresponding to the previously defined 
field ellipsoid. Figures0]|5]and|6]illustrate the density distribu- 
tion of Galactic space velocities projected on the (V, U), (W, V) 
and (W, U) planes, respectively. The uncorrected distributions 
are shown in the upper panels, while the subtracted Gaussian 
distributions and the corrected density distributions are pre- 
sented in the lower panels. 
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Fig. 5. Density distribution of Galactic space velocities pro- 
jected on the [W, V) plane for the whole comparison sample. 
The uncorrected density distribution, the subtracted Gaussian 
distribution and the corrected density distribution are presented 
in the upper, central and lower panels, respectively. The field 
ellipsoid is overplotted with the blue line, while the Hyades, 
Sirius and Hercules ellipsoids are indicated by the red, blue 
green and yellow lines. 

The next step consisted in defining the ellipsoids describ- 
ing the remnant "structures" in the corrected density distribu- 
tion. As a first aproximation, we adopted the kinematic param- 
eters of the Hyades, Sirius and Hercules streams derived by 
Fama e~et alJ (2005a), and then adjusted by eye the ellipsoid 
centroids and semiaxes according to what we observed. Tabled 
lists the kinematic parameters finally adopted for each veloc- 
ity ellipsoid. In Figures 0] 13 and |6] the projected ellipsoids are 
represented and identified as the Hyades, Sirius and Hercules 
streams. The projected field ellipsoid is also overplotted in each 
figure (blue line). The number of targets falling into each ve- 
locity ellipsoid are listed in Tabled 
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Fig. 6. Density distribution of Galactic space velocities pro- 
jected on the (W, U) plane for the whole comparison sample. 
The uncorrected density distribution, the subtracted Gaussian 
distribution and the corrected density distribution are presented 
in the upper, central and lower panels, respectively. The field 
ellipsoid is overplotted with the blue line, while the Hyades, 
Sirius and Hercules ellipsoids are indicated by the red, blue 
green and yellow lines. 



The method used to obtain the field ellipsoid does not make 
any a priori assumptions concerning the dynamical streams. A 
different approach would be to adopt a negligible vertex devi- 
ation for the field ellispoid by conside ring that the Hercules 
stream is its main contributor (e.g. iMiihlbauer & Dehnenl 
2003). We have checked that the assumption of a negligible 
vertex deviation does not change our results. Although the 
Hercules ellipsoid results slightly enlarged, and consequently 
the number of comparison stars belonging to the stream is en- 
hanced from 5 to 8, no planet-host stars are found within the 
modified Hercules ellipsoid. 
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Next, we defined the groups in which our sample would 
be divided for the comparison: planet hosts, comparison sam- 
ple, metal-rich comparison stars and comparison stars with so- 
lar metallicity or lower. The metallicity cutoff was chosen so 
that the median of the metal-rich comparison sample was equal 
to the one of the planet-host distribution ([Fe/H]=0.15, see 
Figure |5J upper panel), and also so that the mass-metallicity, 
or age-metallicity, relation affected as few targets as possible 
(see Figure|5J lower panel). Cutting the comparison sample at 
[Fe/H] = 0.07 fulfilled the condition concerning medians. The 
metallicity cutoff was then rounded off to [Fe/H] =0.1. 

We then proceeded to count how many stars from each 
group fell into the ellipsoids corresponding to the field and 
to the Hyades, Sirius and Hercules streams. We also counted 
the numbers expected in the streams if the group were to have 
followed a Gaussian distribution, with <x equal to the veloc- 
ity dispersion of the group. For the planet-host sample, we did 
this for two Gaussian distributions, the first with the dispersion 
characteristic of the whole comparison sample, and the second 
with the dispersion characteristic of the metal-rich comparison 
subsample. The count for each stream was corrected by sub- 
tracting the expected number from the Gaussian distribution 
and normalized by dividing by the total number of the group in 
the field. This corresponds to the ratio of excess objects in the 
stream. The Sirius and Hercules streams were excluded from 
this analysis since the two velocity ellipsoids contained respec- 
tively one and no planet host stars. 

A further analysis was carried out by means of the density 
distributions of the Galactic space velocities, derived by using 
the fifth nearest-neighbour method. For each group indepen- 
dently, we computed the density distribution corresponding to 
a Gaussian distribution with <x equal to the velocity dispersion 
of the group. In the case of the planet-host group, two Gaussian 
distributions were tested, one with the dispersion of the whole 
comparison sample, and the other with the dispersion of the 
metal-rich comparison subsample. The "observed" distribution 
was then corrected by the Gaussian contribution. From the cor- 
rected density distribution, we derived an average density for 
the field and the Hyades ellipsoids. The averaged field density 
was then subtracted that of the Hyades in order to obtain the 
overdensity over the mean field "level". 

3. Results and discussion 

3.1. Thin/thick disc 

Our results reveal that the CORALIE sample bears one planet- 
host target, HD4308, and 21 stars with no known planetary 
companions, with kinematics typical of the thick disc (see 
FigureG}. The metallicities of HD 4308 ([Fe/H] = -0.27 from 
the CCF calibration) and of the thick disc targets without 
known planets are typical of the thick disc population (see 
Figure^). Only one comparison star, HD 152391, shows a su- 
persolar metallicity ([Fe/H] = 0.03 from the CCF calibration). 
The following study combining both kinematical and chemi- 
cal information will be limited to the subset of CORALIE stars 
with available ho mo geneous abundanc es from lEcuvillon et ail 
(2004a b, 2006b) and lGilli et al1 (E()06). 
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Fig. 7. Toomre diagram for the CORALIE targe t s with c hem- 
ical abundances by lEcuvillon et alJ J2004albL l2006bl) and 
iGilli et alJ (H)06). Dotted lines indicate constant peculiar ve- 
locities v pec = (t/^ SR + V^ SR + Wlsr)' /2 in steps of 50kms _1 . 
Filled and open symbols represent stars with and without plan- 
ets, respectively. Targets that are full and intermediate members 
of the thick disc population are indicated with red triangles and 
squares, respectively, while full and intermediate members of 
the thin disc are represented by blue diamonds and squares, re- 
spectively. The four comparison targets full members, and the 
four planet-host stars intermediate members of the thick disc, 
are identified by their numbers. 



The Toomre diagram for the subset of C ORALIE stars 
with c hemi cal abunda n ces fr om lEcuvillon et alJ ( 2004a b , 
l2006bl) and lGilli et alJ J2006I) is shown in Figure Five 
comparison stars are very likely members of the thick disc 
population (HD 20794, HD 135204, HD 136352, HD 152391, 
HD 191408), while no planet-host targets present such a high 
probability. Only the four planet hosts HD6434, HD 13445, 
HD 1 1 1232 and HD 1 14729, are potential members of the thick 
disc, with probabilities P(thick disc)/P(thin disc) between 1 
and 10. The [X/Fe] vs. [Fe/H] trends for oxygen and magne- 
sium are shown in Figure|S] and in Figures|5]and[TOlfor the rest 
of the elements (C, N, Na, Al, Si, S, Ca, Sc, Ti, V, Cr, Mn, Co, 
Ni, Cu and Zn). 

Previous works (e.g. Fuhrmann 1998: iFeltzing et alll2003l: 
iBensbv et all2003l) reported different abundance trends for the 
thin and thick disc populations. In particular, the cr-elements 
oxygen and magnesium show clearly distin ct trends for the thin 
and thick discs at sub-solar metallicities jBensbv et alJl2003l 
2004). However, IChen et alJ d2000l) presented results from a 
sample of 90 F and G dwarfs which show no significant scatter 
in a-element ratios as function of [Fe/H]. Another interesting 
issue is the existence of stars with thick disc kine matics at high 
metallicities, point ed out by several authors (e.g. lBensbv et alJ 
l2003UMishenina et al.l2004l) . 

Our results confirm the possibility of thick disc stars at 
[Fe/H] > 0, since one of the five comparison stars with 
thick disc kinematics that was previously chemically analysed, 
HD 152391, has [Fe/H] = 0.03. The same target has been al- 
ready reported as thick disc member with supersolar metallic- 
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Fig. 8. rO/Fel and [M g/Fel vs. [Fe/H] for CORALIE targets with chemical abundances bv lEcuvillon et all J2004albl l2006b) 
and lGilli et all d2006l) . Filled and open symbols represent stars with and without planets, respectively. Targets which are full 
and intermediate members of the thick disc population are indicated with red triangles and squares, respectively, while full and 
intermediate members of the thin disc are represented by blue diamonds and squares, respectively. 



ity in the kinematic study of the whole CORALIE sample (see 
the first paragraph of this section). It is interesting to note how 
the metallicity derived by the CCF calibration agrees perfectly 
with the iron abundance iss ued by the detailed spectroscopic 
study bv lSantos etaD J2005t) . 

Nevertheless, we did not find any difference in the abun- 
dance patterns marked by thick and thin disc stars, nor any 
significant signature of enrichment in a-elements for thick disc 
population. For magnesium (see Figure|8j right panel), it might 
be possible that the thick disc members lay on the upper en- 
velope of the points. However, the small number of possible 
members for this population in our sample does not enable us 
to reach any conclusion. 

3.2. Dynamical streams 

Table[3]shows the results obtained with the "counting" method, 
using the ellipsoid definitio ns from our work (Definition 7), 
and the ones from iFamaev et all J2005bh (Definition 2). The 
first column lists the different groups: planet-host stars, assum- 
ing for the Gaussian the cr of the whole comparison sample, 
and of the metal-rich subsample; the whole comparison sam- 
ple; and the metal-rich and the solar-metallicity comparison 
subsamples. Then are listed the numbers of targets in the field 
and Hyades ellipsoids, the expected number of targets from the 
Gaussian distribution, and finally the ratio of objects in excess 
from what expected for a Gaussian distribution. The uncertain- 
ties were estimated as Poissonian errors. 

The results issued from this method indicate that the 
dumpiness typical of the Hyades is much more significant 
for the metal-rich subsample than for the whole sample. This 
confirms the high average metallicity report ed for Hyades 
stream members by prev ious studies (e.g. lRaboud et alJll998t 
IChereul & Grenonl200ll) and agrees with the mechanism pro- 
posed as an explanation of the anom alies in the (U, V) plane 
JPehnenlEoOot IFamaev et alJl2005al) . According to this de- 
scription, the stars would have formed at inner Galactic radii, 
with higher metallicity, and then be trapped into resonant or- 



bits by the non-axysimmetric potential, and thrown closer to 
the Sun. 

The planet-host group shows an excess in the Hyades 
stream much more significant than the other groups. Errors as- 
sociated with stars with planets are large because of the small 
number of targets with planetary companions involved in the 
whole volume-limited sample CORALIE. However, the targets 
with planets present an excess in the Hyades stream signifi- 
cantly larger than that obtained for the other groups. In particu- 
lar, the dumpiness of planet-host stars is different from that of 
the whole comparison sample, with a significance of more than 

1.5 cr. 

The kinematic behaviour of CORALIE stars with planets 
seems to match better with that observed in the metal-rich com- 
parison subsample than with the whole comparison sample. 
Both metal-rich and planet-host samples show a much larger 
dumpiness than the whole comparison sample. Nevertheless, 
it is important to note that the presence of planet-host stars in 
the Hyades stream is even denser than that found for metal- 
rich "single" stars. This outcome is not altered by changing 
the assumption of the characteristic velocity dispersion of the 
group with planets; there is no significant difference in adopt- 
ing a Gaussian with the cr of the whole comparison sample, or 
of the metal-rich subsample. 

If the kinematic parameters of dynamical streams by 
IFamaev et alJ J2005al) are adopted (see Table [3] Definition 2), 
the results are similar to those reported above. The differences 
between the excesses for the comparison groups obtained us- 
ing the two kinematic parameter sets are mostly irrelevant. 
For planet-host stars, the two definitions give larger discrep- 
ancies. In the Definition 2 case, the group of stars with plan- 
ets behaves in the same way as the metal-rich comparison 
subsampl e. However , it is worth noting that the analysis by 
Fam aev et alJ d2005ah involved K and M giants, while the tar- 
gets of the volume-limited CORALIE survey are mostly G and 
K dwarfs. In that sense, the determination of the dynamical 
stream parameters could be affected by the target spectral type. 
The kinematic parameters of the dynamical streams obtained 
in our work are thus considered more suitable for our analysis. 
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and Gilli et al. (2006). Filled and open symbols represent stars with and without planets, respectively. Targets which are full 
and intermediate members of the thick disc population are indicated with red triangles and squares, respectively, while full and 
intermediate members of the thin disc are represented by blue diamonds and squares, respectively. 



The results from the "overdensity" method are listed in field density to the Hyades overdensity. In the same table we 

Table|4] The Def. 1 column refers to the case adopting the kine- also included the same results for different cases: 
matic parameters of the streams ellipsoids from our work and 

[Fe/H] = 0. 1 as metallicity cutoff, and subtracting the averaged n c ■ . • , • . • . t 

L ' 1 J ' & fe - Def. 2: adopting the kinematics parameters of the stream 

ellipsoids bv lFamaev et alJ J2005b): 
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- Def. 1 mod. 1: adopting our kinematic parameters, but with- 
out subtracting the averaged field density; 

- Def. 1 mod. 2: adopting our kinematic parameters, but 
adopting [Fe/H] = 0.05 as metallicity cutoff between the 
metal-rich and solar comparison subsamples; 

- Def. 1 M < 1.2M : only the targets with M < 1.2M were 
included in the analysis in order to check the effect of age. 



The errors have been estimated by the bootstrap method, where 
the estimate of standard error is the standard deviation of the 
bootstrap replications. 

The excess values obtained from the "counting" method 
(see Table|5J Definition 1), and the overdensities obtained from 
the "overdensity" analysis (see Tabled Defl), and their cor- 
responding uncertainties, are represented in Figure ^2 1° of- 
fer a more immediate representation. The black and red trian- 
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Table 3. Results for the "counting" method in the Hyades stre am. The Defin i tion 1 a nd Definition 2 columns refer to the adopted 
definitions for the stream ellipsoids, from our work and from lFamaev et alJ (12005a). respectively. The numbers of stars lying in 
the field and Hyades ellispoids are listed, as well as the expected number of stars for a Gaussian distribution, and the resulting 
ratio of objects in excess in the Hyades stream. The assumed uncertainties are Poissonian errors. 



Group description 


Field 


Hyades 


Definition 1 
Exp Hyades 


Excess 


Field 


Hyades 


Definition 2 
Exp Hyades 


Excess 


Planet-hosts: cr whole comp. 


47 


10 


1.7 


0.18±0.07 


47 


7 


1.4 


0.12+0.05 


cr M-rich comp. 






2.3 


0.16±0.06 






1.9 


0.11+0.06 


Whole comparison 


440 


40 


16.3 


0.05+0.01 


440 


33 


13.0 


0.05+0.01 


M-Rich comparison 


119 


19 


5.7 


0.11 ±0.04 


119 


17 


4.7 


0.10+0.03 


Solar comparison 


321 


21 


11.9 


0.03+0.02 


321 


16 


9.2 


0.02+0.01 



Table 4. Results for the "overdensity" method in the Hyade s stream. The Def.l and Def.2 columns refer to the adopted definitions 
for the velocity ellipsoids, from our work and from Fa maev et alJll2005al) . respectively. Def.l mod.l corresponds to the case with- 
out substracting the average density of the field, while Def.l mod.2 reports the results for a metallicity cutoff [Fe/H] c „ /O y/=0.05 
(see the Subsection l3.2l for more details). The last column, Def.l M<1.2M G , refers to the results obtained by excluding the targets 
with M>1.2M . The uncertainties were derived applying the bootstrap method. 









Hyades Overdensity 




Group description 


Def. 1 


Def. 2 


Def. 1 mod. 1 


Def. 1 mod. 2 


Def. 1 M < 1.2M G 


Planet-hosts: cr whole comp. 


0.031+0.017 


0.030+0.017 


0.031 


0.031+0.018 


0.011+0.020 


cr M-rich comp. 


0.029+0.018 


0.027+0.018 


0.028 


0.028+0.017 


0.009+0.020 


Whole comparison 


0.004+0.004 


0.004+0.004 


0.004 


0.004+0.004 


0.003+0.004 


M-Rich comparison 


0.020+0.011 


0.019+0.012 


0.020 


0.012+0.007 


0.013+0.011 


Solar comparison 


0.003+0.005 


0.003+0.004 


0.003 


0.002+0.003 


0.003+0.005 



gles represent the results for planet-host stars, corrected with a 
Gaussian distribution with cr equal to the velocity dispersion of 
the whole comparison sample and of the metal-rich subsample, 
respectively. 

The "overdensity" analysis (see Table|4] Def.l) leads to re- 
sults similar to the "counting" method (see Figure ITTb . On the 
one hand, the comparison metal-rich subsample shows a more 
significant overdensity in Hyades with respect to the whole 
comparison sample, agreeing with what reported by previous 
studies. On the other hand, the behaviour of planet-host stars is 
consistent with what we observe for the metal-rich CORALIE 
"single" stars, whereas it does not seem to reproduce the kine- 
matic characteristic of the whole comparison sample. In partic- 
ular, stars with planets present an overdensity in Hyades even 
more evident than the one found in the metal-rich comparison 
subsample. Planet-host overdensity stands out from the value 
for the comparison sample in 1 .5 cr. 

No differences em erge if we adopt the stream parameters by 
iFamaev et alJ j2005a) (see TableH Def.2), or if we do not sub- 
tract the average density of the field (see Table|4] Def.l mod.l). 
For the case adopting [Fe/H] = 0.05 as metallicity cutoff (see 
Table |4] Def. 1 mod. 2), there are no significant changes in the 
results, except for the overdensity of the metal-rich compari- 
son subsample, which slightly diminishes. In all the cases, the 
planet-host stars are much denser in the Hyades stream with 
respect to the whole comparison sample, with a significance of 

1.5cr. 



If younger stars are excluded by considering only targets 
with M < 1.2 M Q , the results for the group of stars with 
planets deviate from the previous ones (see Table |4] Def. 1 
M < 1.2 M Q ). In this case, planet-host stars show a kine- 
matic behaviour almost indistinguishable from the metal-rich 
comparison subsample, with an equivalent overdensity in the 
Hyades. Younger planet-host stars might be thus responsible 
for the outstanding excess of stars with planets with respect to 
the rest. 

This raises a question concerning the possible effects of age 
on the kinematic conduct of the different groups studied. The 
small number of targets with M > 1 .2 M Q does not allow us 
to study the behaviour of the younger tail independently, so no 
conclusive argument can be added to discard this possibility. 
All we can say is that, roughly, a higher percentage (around 
50%) of younger stars is found among planet-host targets with 
Hyades stream kinematics than among comparison stars of the 
Hyades stream. 

However, the ratio of younger stars in the planet-host group 
is the same as that found in the metal-rich comparison subsam- 
ple: 29 younger over 1 15 older metal-rich comparison targets, 
and 1 1 younger over 50 older planet-host stars. This ascertain- 
ment permits us to be more confident that the results are not 
due to a larger percentage of younger dwarfs in the planet-host 
sample than in the other groups. The age effect for planet-host 
targets seems to be strictly related to their membership to the 
Hyades stream. Moreover, this seems to confirm that planet- 
host stars behave similarly as the metal-rich population of the 
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Fig. 11. Results from the "counting" (left panel) and the "overdensity" methods [right panel). Black, red and blue diamonds 
represent the whole comparison sample, and the metal-rich and solar subsamples, respectively. Black and red triangles indicate 
planet-host stars, adopting a Gaussian distribution with <x equal to the velocity dispersion of the whole comparison sample, and 
of the metal-rich subsample, respectively. The error bars show Poissonian errors for the "counting" method (left panel), and 
uncertainties derived by the bootstrap method for the "overdensity" analysis (right panel). 



solar neighbourhood. It thus represents a further argument in 
favour of a primordial origin of the metal excess found in stars 
with giant planets. It cannot be discarded that the important 
dumpiness of young planet-host stars in the Hyades stream 
might be due to a more efficient planet formation in the pri- 
mordial Hyades cluster, progenitor of a significant fraction of 
young members of the Hyades stream. 

All our results suggest that planet host stars follow the 
kinematic conduct of the metal-rich comparison subsample, 
showing an overdensity in the Hyades stream with respect to 
the whole comparison group. Furthermore, the younger planet- 
host targets might present a more significant dumpiness in the 
Hyades stream that would result in a further overdensity in the 
Hyades stream, exceeding that found in the metal-rich com- 
parison subsample. However, much more targets are needed to 
confirm this latter suggestion. 

This supports the "primordial" hypothesis, according to 
which the metal-rich nature of planet-host stars would be due to 
the high metal content of the protoplanetary clouds from which 
the systems formed. Moreover, it can also be understood in the 
scen ario proposed as a po ssible origin of the dynamical streams 
(see lFamaev et all2005bl and references therein). According to 
this hypothesis, the kinematic members of the Hyades stream 
would have formed in Galactic inner regions, where the inter- 
galactic medium is more metal-rich. Their orbits would have 
been successively perturbed by spiral waves, causing radial mi- 
gration towards the solar neighbourhood. 

The fact that ljames et alJj200 6) did not find any metal-rich 
population in a sample of three nearby young stellar forming 
regions could support the scenario in which most metal-rich 
field stars come from inner regions in the Galaxy. Moreover, 



Chiappini in preparation) found that Galactic chemical 

evolutionary models, which reproduce most of the abundance 
trend peculiarities for thin and thick disc populations, cannot 
explain the existence of the most metal-rich stars of the solar 
vicinity (at ~ 8 kpc). However, those objects are consistent if 
the same models are computed at inner Galactic radii (~ 4 kpc). 

iDehnenl |2000) proposed the Outer Lindblad Resonance 
(OLR) as the mechanism reproducing the bimodality of the lo- 
cal velocity distribution, identified as the Hercules stream. In 
this model, closed orbits inside the OLR would move slightly 
outwards, while those outside the OLR would move inwards, 
producing two stellar streams. Since the simulations show that 
the radius of the OLR of the Galactic bar lies inside the so- 
lar circle, the stellar velocity distribution observed for late-type 
stars in the solar neighborhood would be affected by the OLR. 
Thes e results were then con firmed by further simulations by 
Miihlbauer & Dehnen This has been the first example 

of a non-axysimmetric origin for a stream and has opened up 
new possibilities for relating the other streams with other non- 
axysimmetric perturbations. 

Several simulations have studied the effects of non- 
axysimmetric perturbations, in particul ar of transient spi 
ral waves, on the galactic disc (e.g. ISellwood & Binnev 



|2002 [lBissantz et alJ2003t lQialfl-aba rt"vi2004lDe Simone et al 
2004). They found that spiral waves can cause radial migration 
near the corotation radius. According to Sellwood & Binney's 
model, stars just inside corotation would swap places with 
those outside, and bot h groups just would exchange places. 
In this scenario, Famae vet alJ J20 05a) interpreted Hyades and 
Sirius streams as the inward- and outward-moving streams of 
stars on horshoe orbits that cross rotation. This would eas- 
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Fig. 12. Density distributions of Galactic space velocities projected on the (V, U) plane, for the metal-rich (left panel) and solar 
(right panel) comparison subsamples. The uncorrected density distribution, the substracted Gaussian distribution and the cor- 
rected density distribution are presented in the upper, central and lower panels, respectively. The field ellipsoid is overplotted 
with the blue line, while the Hyades, Sirius and Hercules ellipsoids are indicated by the red, blue green and yellow lines. 



ily explain the seemingly peculiar metallicity of these moving 
groups. 

Another interesting feature appears if we inspect and com- 
pare the density distributions of Galactic space velocities for 
the metal-rich and solar subsamples (see Figure I12> . In the 
Hyades ellipsoid, the peaks corresponding to the overdensities 
of the two groups are displaced: the overdensity of the metal- 
rich subsample (in the left panel) is shifted to lower U val- 
ues, whereas the solar subsample (right panel) shows a peak 
shifted towards higher U velocities. It is also worth noting that 
the Sirius overdensity completely disappears when considering 
only metal-rich targets (see right panel). This suggests that the 
members of the Sirius stream have solar and sub-solar metallic- 
ities. The refore, it is in good agreement with the scenario pro- 
posed by Fam aev et alJ ( 12005 al) . according to which the Sirius 



stream would correspond to the inward-moving counterpart of 
the Hyades stream. 

In this framework, the larger overdensity of planet-host 
stars observed in the Hyades stream would suggest that at in- 
ner galactic radii, in environments more metal-rich than the so- 
lar neighborhood, systems with giant planets form more easily. 
This strongly supports the primordial origin proposed to ex- 
plain the metal e xcess o bserved in stars harbouring giant plan- 
ets JSantos etalJl2000l 1200 ll) . Current models of planet for- 
mation by core accretion are in agreement with this observa- 
tional constraint, and find that m etal-rich systems favour the 
forma tion of massive planets (e.g. llda & LinlEo05l iBenz et alJ 
2006). The combination of this effect with the detection bias 
towards more massive planets enables models to reproduce the 
observed correlation betwee n stellar m etallicity and the likeli- 
hood of harbouring planets tBenz et all 2006). 
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4. Concluding remarks 

This is the first study carried out on a volume-limited sam- 
ple of stars with and without exoplanets aimed at investigating 
the comparative kinematics of planet-host stars and their pos- 
sible relation with the dynamical streams of the solar neigh- 
bourhood. With this goal, we have compared the group of stars 
harbouring planets with the rest of the CORALIE sample, dis- 
tinguishing subgroups of metallicity and age. This was made 
possible by the precise radial velocity measurements and the 
CCF parameters included in the CORALIE database, and using 
the HIPPARCOS and Tycho-2 catalogues. For targets with de- 
tailed abundance determinations, a complete study on whether 
they belong to the thick and thin disc populations, and on pos- 
sible signatures in the abundance trends, has been carried out. 
The analyses related to the dynamical streams have been car- 
ried out using two independent methods that have produced 
equivalent results. 

The main outcomes of this work are the following: 

- The CORALIE planet-host target HD 4308 has been found 
to be a very probable member of the Galactic thick disc, to- 
gether with another 21 comparison stars; the dependence of 
the abundance trends on the thin/thick disc populations has 
been studied for a subset of the CORALIE sample with ho- 
mogeneous abundance determinations from previous stud- 
ies. 

- No significant differences or peculiar enrichment in a- 
elements are observed in the abundance patterns of the sub- 
set of stars with homogeneous abundance determinations 
from previous studies. However, the small number of thick 
disc members prevents us from reaching conclusive results. 

- The tentative existence of thick disc stars at supersolar 
metallicities is confirmed by the appearance of one thick 
member over 21 with high metallicity (HD 152391 with 
[Fe/H] = 0.03). 

- The high average metallicity of the Hyades stream reported 
by previous works is confirmed. Our results are in agree- 
ment with the scenario of non-axysimmetric perturbations, 
such as transient spiral waves, proposed by several results 
of simulations and observations. According to this model, 
the Hyades stream would be the outward-moving stream of 
high-metallicity stars born at inner Galactic radii and then 
pushed into the solar neighborhood. We have also found a 
new clue confirming the solar and subsolar average metal- 
licity suggested for the Sirius stream. 

- The group of planet-host stars shows a kinematic behaviour 
much more similar to the metal-rich comparison subsam- 
ple, than to the whole comparison sample, in the sense that 
the overdensity in the Hyades stream observed for stars 
with planets is much higher than that observed for the 
whole comparison group. This issue has been reached by 
the two methods independently and strongly supports a pri- 
mordial origin for the metal excess observed in stars with 
giant planets. 

- If we interpret our results within the scenario proposed to 
explain the presence and origin of dynamical streams in the 
solar neighbourhood, stars with giant planets could have 
formed more easily at inner Galactic radii in a more metal- 



rich intergalactic medium and then suffered radial displace- 
ments due to a non-axysimmetric component of the galactic 
potential, pushing them into the solar neighbourhood. 

This paper presents the results of the first complete analysis 
on this topic. However, the uncertainties affecting our results 
are large, due to the intrinsic errors of the kinematics analysis, 
but also, and mainly, to the inevitably small number of targets 
with planets in the CORALIE sample. Consequently, substan- 
tial improvements will be possible in the near future, when the 
steadily increasing number of known stars with planetary com- 
panions will have risen enough to permit more definite conclu- 
sions. 
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